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ABSTRACT 

In contrast to the Earth, where frictional heating is typicahy neghgible, we show 
that drag mechanisms could act as an important heat source in the strongly-forced 
atmospheres of some exoplanets, with the potential to alter the circulation. We modify 
the standard formalism of the atmospheric energy cycle to explicitly track the loss of 
kinetic energy and the associated frictional (re)heating, for application to exoplanets 
such as the asymmetrically heated "hot Jupiters" and gas giants on highly eccentric 
orbits. We establish that an understanding of the dominant drag mechanisms and 
their dependence on local atmospheric conditions is critical for accurate modeling, not 
just in their ability to limit wind speeds, but also because they could possibly change 
the energetics of the circulation enough to alter the nature of the flow. We discuss 
possible sources of drag and estimate the strength necessary to significantly influence the 
atmospheric energetics. As we show, the frictional heating depends on the magnitude 
of kinetic energy dissipation as well as its spatial variation, so that the more localized 
a drag mechanism is, the weaker it can be and still affect the circulation. We also use 
the derived formalism to estimate the rate of numerical loss of kinetic energy in a few 
previously published hot Jupiter models with and without magnetic drag and find it to 
be surprisingly large, at 5-10% of the incident stellar irradiation. 

1. Introduction 

While the field of atmospheric dynamics is well established and mature, we are currently 
watching its new expansion into the exotic regimes introduced by extrasolar planets. Much of 
the same, decades-old theory can be fruitfully applied to exoplanet atmospheres, but important 
differences exist between these planets and those found in our solar system. This may result in the 
need to adapt theories to a new context and check that basic assumptions are not applied where 
they are not valid. 
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The strongly-irradiated atmospheres of close-in, extrasolar gas giants ("hot Jupiters") differ 
from solar system giants in a few important ways: 1) the stellar heating of the atmosphere dominates 
over internal heating, 2) these planets are expected to be tidally locked into synchronous orbits, 
meaning that the irradiation pattern is strongly asymmetric, and 3) the long rotational (= orbital) 
period means that the Coriolis force will have a weaker effect in these atmospheres (compared to 
Jupiter or Saturn) and atmospheric features should exist on a planetary scale (for a comprehensive 
review see Showman et al. 2010). One important issue that is currently being investigated is the 
role of friction or drag in these atmospheres. 

Atmospheres in the solar syst em are largely iriyiscid and estimates for hot Jupiters also imply 



very high Reynolds numbers (e.g. iLi &: Goodmanll20ld ): h owever, there can be other important 



sources of friction or drag in an atmosphere. In particular, iGoodmanI ()2009l ) questioned whether 
sources of drag were being correctly accounted for in hot Jupiter numerical models, concerned that a 
lack of drag could allow the winds to accelerate to unphysical speeds. On the Earth it is well known 
that surface drag has an important impact on the momentum budget of the atmospheric flow. While 
hot Jupiter atmospheres will not experience drag from a solid surface, other relevant mechanisms 
may be at work (e.g., shockj^, instabilities, and magnetic drag, which we review in more detail in 
Section [3.ip . It may be that models without explicit friction are still valid, if physical sources of 
friction are weak (compared to the radiative forcing) and/or numerical dissipation dominates, but 
it is necessary to have good estimates of the strengths of these mechanisms to determine if that is 
the case. 

In this paper we will discuss the importance of friction in the context of the energetics of 
atmospheric circulation. We emphasize that drag can have multiple effects on the flow, namely as: 
1) a momentum sink, or 2) as a possible source of the heating that drives the flow. (Throughout 
this paper we will be using the terms friction and drag synonymously. Whether it is frictional 
dissipation from microscopic interactions, or large-scale drag from other physical processes, the 
result will still be the irreversible conversion of kinetic energy into localized heating.) While the 
concerns raised in lCoodmanl (j2009l ) were related to the first point, in this paper we will show that 
the second point may be of comparable importance in terms of correctly modeling exotic exoplanet 
atmospheres. 

In Section [2] we review standard equations for the global energy balance in an atmosphere 
(from Pearce 1978) and explicitly account for heating from frictional dissipation. We then discuss 
the contexts in which it is acceptable to neglect friction from the energetics (even if not from the 
momentum equation). In Section [3] we discuss the difference between these contexts and the case of 
a strongly-forced atmosphere, with a discussion of the drag mechanisms that have been proposed 
for hot Jupiters (13. ip . We then present equations for calculating the local values of the various 
energies and their conversion rates in Section 13. 2( and in Section 13.31 use these to estimate whether 
frictional heating should have been explicitly accounted for in previous numerical models. We 



^We are loosely referring to shocks as drag, in the sense that they transform bulk kinetic energy into heat. 
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briefly discuss how these concepts will apply to the special case of an eccentric, close-in gas giant 
planet in Section HI Finally, in Section [5] we summarize our main points and emphasize the context 
in which friction may matter to the energetics of atmospheric circulation. 



2. The Standard Theory of Friction in Atmospheric Energetics 



The total energy in an atmosphere is divided between several forms: the specific kinetic energy 
where v is the three-dimensional velocity), the internal energ}!^! (c^T, where is the specific 
heat and T is the local temperature), and the gravitational potential energy [gz, where g is the 
gravitational acceleration and z is altitude). If the atmosphere is in hydrostatic equilibriurcQ 
{dP/dz = —pg, where P is pressure and p is density), the ideal gas law (P = pRT) and the 
relation between the specific gas constant and specific heats {R = Cp — c^) can be used to show 
that the integrated internal and gravitational potential energies for a column of air are equal to the 
integrated total specific enthalpy (cpT). In studying the energetics of an atmosphere, we want to 
track the sources and sinks of potential energj|§ (through heating and cooling), the kinetic energy 
loss (through frictional dissipation), and the adiabatic conversion between potential and kinetic 
energies. This conversion occurs when parcels of fluid flow down gradients created by differential 
heating in the atmosphere. 



The concept of "available" potential energy (APE) was first formulated by iLorenj (jl955l ) as 
the component of energy in an atmosphere that is generated by heating and available for direct 
conversion into kinetic energy. His definition was in refer e nce t o some idealized, APE=0, state 
with no horizontal temperature/pressure gradients. iPearcd (| 19781 ) created a new formalism based 
on sources, sinks, and conversions, in which he recognized that APE is only generated through 
variations in heating, whereas uniform heating only creates "unavailable" potential energy (UPE) . 
This coricept of APE has also been form ulated a s avai lable exergy or enthalpy (jPuttonl Il973l : 
Marquet][l99ll ). We find the formalism by iPearcd (|l978l ) most useful for our purposes here and 
in this section will follow his derivations. Given that these are generally unfamiliar concepts in 
astrophysics, we consider it useful to review them here. 



If X is a local quantity (per unit mass) in the atmosphere, its global average is denoted by X 



''The internal energy of a parcel of air can change through compositional or phase changes (e.g. evaporation or 
condensation), but in this analysis we assume a constant state. 

^This is a common assumption — and usu ally justifiably so— for most atmospheres. For a discussion of hydrostatic 



equilibrium in hot Jupiter atmospheres see IShowman et al.l (|2008l ). Throughout this work we will assume that the 



atmospheres we consider are in hydrostatic equilibrium. 

^Note that throughout this paper we will be using the term "potential energy" to refer to thermodynamic potential 
energy (usually in the form of enthalpy). For an atmosphere in hydrostatic equilibrium the gravitational acceleration 
is balanced by the vertical pressure gradient; it is not this component of the potential energy that leads to winds, 
but rather the horizontal gradients of potential energy. 
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and calculated as: 

^" = if / = i? / f ""^"^ = ^ / ^"""'"'^ 

where is the total mass of the atmosphere, x and y are the horizontal coordinates, and the vertical 
coordinate has been converted from altitude to pressure by assuming hydrostatic equilibrium and 
that the gravitational acceleration (g) is constant. 



Following iPearcd (jl978l ) , we can then write global expressions for the time derivatives of the 



specific total potential (etot) and kinetic (bk) energies of the atmosphere: 

detot/dt = -c + q (2) 

and 

dex/dt = c — d (3) 

where q is the globally averaged heating rate, c is the rate of conversion from potential to kinetic 
energy, and d is the rate of dissipation of kinetic energy. (In principle the conversion could go in 
either direction, but the loss of kinetic energy through dissipation requires a positive net conversion 
from potential to kinetic energy.) To evaluate what fraction of the total potential energy will be 
available for conversion into kinetic energy, we begin with the thermodynamic equation: 

ds/dt = {cp/T){dT/dt) - Ruj/P = q/T (4) 

where s is the specific entropy and uj = dP/dt. The global average of this equation is: 

dsM = (I) = g (1) + <Z* (5) 

where * denotes the deviation from the global average. We can then define a reference temperature, 
Tr, by 1/Tr = and multiply Equation [5] by Tr to obtain: 

Trds/dt = q + q* = du/dt (6) 

The quantity T^ds/dt is, by definition, the change in energy associated with a globally uniform 
change in entropy and therefore defines the time derivative of the specific "unavailable" potential 
energy (UPE, or u). The potential energy that is available for conversion to kinetic energy (the 
APE, or a) will be the difference between this and the entire potential energy of the atmosphere 
(from Equation ^ : 

dd/dt = detot/dt - du/dt = -c- q*(Tyf)* (7) 



Through this derivation iPearcd (|l978l ) was able to identify the term —q*(Tr/T)* (the global average 



encompasses all variables) as the efficiency with which heating acts as a source for APE rather than 
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UPE, and see that it depends on the spatial correlation between non-uniform heating and a non- 
uniform temperature structurejf] 

We now return to the equations for the total potential and kinetic energies ([21 [3]) to see that 
a steady state atmosphere {detot/dt = dex/dt = 0) requires: 

q = c = d (8) 

In order to explicitly track the frictional heating due to the dissipation of kinetic energy, we separate 
the heating rate into frictional and non- frictionaE components, Q = Qf + Qrad- We can then equate 
the global loss of kinetic energy with the global frictional heating (d = qf , an irreversible conversion) 
and see that Equation [8] implies: 

q = Qf, Qrad = 0. (9) 

The expression grad = is simply a statement of global radiative balancelfl We can also write the 
heating deviation as a combination of radiative and frictional components: 

Q* = Qrad + Q*f- (10) 



Taken all together, these equations from iPearcd (|l978l l identify the sources of potential energy in 



a steady-state atmosphere in global radiative equilibrium: 

• The only source of UPE is the mean frictional heating. 

• The source of APE is differential heating, both radiative and frictional. 



With these equations we can now clearly formulate the context in which frictional heating 
can be neglected from the atmospheric energetics. If the friction is mostly uniform it will all go 
into generating unavailable potential energy and will not directly affect the circulation. If there 
is a good estimate of the amount of non-uniform frictional heating, it can be compared to the 



Throughout this analysis we are implicitly assuming time-averaged values, including the definition of the refer- 
ence temperature, Tr, as is appropriate for atmospheres with minimal global-scale time variability. However, these 
equations are also valid for instantaneous values and could be expanded to analyze the energetics of a time-varying 
heating and temperature structure. 

^To be more explicit, we make a distinction between frictional heating — the increase in molecular kinetic energy 
from frictional dissipation of the kinetic energy of the flow — and non-frictional heating, which includes all other sources 
of heating, including radiation, conduction, and latent heat release. Typically the radiative heating dominates over 
other non-frictional sources. 

*By assuming that the total energy of the atmosphere is conserved {detot/dt — 0) we are neglecting the evolutionary 
cooling of the planet. The gas giants in our solar system are not in global radiative balance, with all but Uranus 
radiating more flux than they receive from the Sun. However, here we are considering strongly forced atmospheres, 
where the stellar irradiation dominates over the internal heat flux and the amount by which the planet is out of global 
radiative equilibrium should be negligible. 
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radiative heating to determine whether it can be a significant source of APE. In the case of the 
Earth, for example, it has long been know n that variati ons in the frictional heating are well below 
10% of the variations in radiative heating (jPearcdll978l ). Thus in the Earth-modeling context the 
surface friction is an important momentum sink, but can safely be neglected from the energetics 
without altering the global circulation. We emphasize this point: when the friction is weak, and 
mostly uniform, it is safe to neglect it from the atmospheric energetics without altering the global 
circulation. 



3. Extension to Strongly Forced Exoplanet Atmospheres 

We now turn our attention to extrasolar planets, in particular hot Jupiters. While they should 
not have a solid surface (the main source of friction in a terrestrial atmosphere) , there are several 
drag mechanisms that have been suggested to exist in hot Jupiter atmospheres. Most of these 
possible processes are significantly non-uniform and may be very strong. However, the radiative 
forcing is also intense and a careful comparison of the relative strengths is warranted. 

Before moving on to discuss these atmospheres in more detail, we first point out that there is a 
complication in defining the atmosphere of the planet, as separate from the interior. Unlike terres- 
trial planets, there is no clear boundary at a solid surface, but instead a smoother transition between 
the irradiation-driven (upper) atmosphere, a deeper radiative region that is not directly heated by 
the stellar insolation, and the interior convective zone at pressures of hundreds or thousands of 



bars (e.g., iFortnev &: NettelmannI l20ld ). This makes it more difficult to define a closed system 
and calculate global rates of energy generation and exchange, since there will be some interaction 
with the deeper levels. Although the radiative and advective timescales increase dramatically as 
we descend deeper into the planet, there may be some nonzero net transfer of momentum or heat; 
for example, the downward transport of kinetic energy and subsequent return as an extra internal 



heat flux has been proposed as one mechanism for inflating planetary radii (jShowman &: Guillot 



2OO2I ). These exchange rates may be important for describing the evolution of a planet throughout 
its lifetime, but here we are making the implicit assumption that we are averaging over some time 
period much less than the age of the planet (and that global-scale variability is minimal, consistent 
with many numerical models, including the ones analyzed in Section [3. 3p . 

There is currently no good prescription for how to model interaction between the upper at- 
mosphere and deepest layers of hot Jupiters. Numerical simulations typically assume that the 
modeled atmosphere is a closed box, with energy and momentum conserved. For this analysis we 
will likewise assume that the upper atmosphere is decoupled from the interior-that any exchange 
terms are close to zero and can be neglected. In the case that any exchange with the interior is 
mostly uniform, then the APE generation should not be much effected. However, given the strongly 
day-night asymmetric nature of hot Jupiters, it is worth considering this issue more carefully. Once 
there is a better understanding of how best to treat atmosphere-interior coupling, we will want to 
extend the formalism presented here to include exchange terms. 
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3.1. Discussion of possible drag mechanisms 



Although many drag mechanisms could be at work in a planet's atmosphere, here we re- 
view three main physical processes that have be en identified as possibly producing si gnificant 



drag on hot Jupiters: vertical shear instabilities ( Goodman 



2009 



Li fc Goodman! |2010|'). s hocks 



(jPobbs-Dixon et al.ll20ld : iRauscher Menoull20ld ). and magnetic drag (jPerna et al.ll2010al ). 

While numerical simulations are able to capture large-scale instabilitie^, the computational 
limits on resolution for global circulation models could eliminate any instabilities that are triggered 
on sub- grid scales. Additionally, models that assume hydrostatic equilibrium have no formal ver- 
tical acceleration (vertical velocity exists, but results from the continuity equation) and as such, 
these models are unable to capture vertical shear instabilities, which must be parameterized. As an 
attempt to determine what kind of sub-gri d para meterization needs to be included in hot Jupiter 
global circulation models, iLi GoodmanI (|2010l ) perform a linear analysis and non- linear sim- 
ulations of limited regions of a hot Jupiter atmosphere (two-dimensional, with finite extents in 
longitude and altitude). In their simulations they see weak transient shocks, but the main source 
of kinetic energy dissipation is through recurring Kelvin-Helmholtz instabilities. Although there 
is not yet a clear prescription for how to do so, it will be important to determine how to connect 
processes in localized and global simulations to correctly treat the full range of atmospheric scales 
on hot Jupiters. 

Transonic flow is present in most numerical simulations of hot Jupiter atmospheres. While 
transient shocks m ay exist in the atmosphere, many models contain a shock (or hydraulic iuimj^'^r) 
as a steady feature (tPobbs-Dixon fc LiJ2008l: lshowman et akl boosl. bood : IPobbs-Dixon et al.l[20inl : 
Rauscher Menou 20ld ). This is a localized region of strong dissipation and heating that has a 
markedly different spatial structure than the radiative heating pattern, often appearing on the night 
side o f the planet, a fa c t that may stron gly affect the APE ge neration. The primitive equation mod- 
els of IShowman et al.l (|2008l . bood i and IRauscher fc Menoul (folO) do not include shock-capturing 
schemes and the heating associated with this feature is consistent with adiabatic heating of the 
downward flow that is induced by the strong convergence, rather than from any shock-heating. The 
models in iDobbs-Dixon et al.l (120101 ) do include artificial viscosity and so should be able to cor- 
rectly treat any generation of APE from shocks. In principle the differences between these models 
could demonstrate the impact of standing shocks on the energetics of the circulation, but there are 
other important differences between these modeling approaches and it is difficult to disentangle the 



An example of this i s the horizontal shear instability that is triggered in the shallow hot Jupiter model of 



Menou fc Rauschei] (|2009l ) and which prevents flow speeds from reaching large values 



^"A hydraulic jump is the equivalent of a shock in an incompressible hydr ostatically-balanced flow ; mass and 
momentum fluxes are continuous across the jump, but the energy flux is not (cf. Landau fc Lifshitz 19591 ). In models 
that use the primitive equations of meteorology (instead of the full Navier-Stokes equations), vertical soun d waves 
are fll tere d out a nd ho rizontal ones may or may not be, depending on the boundary conditions imposed. See iKalnav 
( 20021 ) or lVallisI \20od ). 
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physics. 



Finally, recent work implies that the atmospheres of hot Jupiter s may be weakly thermally 



ionized and subjec t to interaction with the planetary magnetic field (jBatvgin fc Stevenson 



201C 



Perna et alJ l2010al ). In particular, the mostly neutral flow may be dragged by the interaction 
between free charges and the magnetic field. This mechanism is strongly affected by geometry — 
both as relates to the orientation of the m agnetic field a nd due to conductivity differences between 



the hot dayside and cold nightside — and iPerna et al.l (|2010al ) estimate that the amount of drag 



could be strong enough to significantly alter the momentum of the flow. The drag is associated 
with an induced magnetic field, itself sustained by electric currents which will dissipate ohmically 
in a non-local way. The result of this complex interaction between the fiuid and the magnetic field 
is that the kinetic energy dissipated will not be locally deposited as heating. Nevertheless, the 
energy may be returned to the flow in a way that alters the energetics of the circulation. A full 
MHD simulation may be required to fully understand how this process can affect the atmospheric 
flow. 

There are also numerical sources of dissipation, although these are usually not well character- 
ized and will vary from code to code, depending on the exact scheme used. A known issue with 
global simulations is the build up of noise on small scales and a common techn ique is to use hyper- 



dissip ation (or hyperdiffusion) to remove power from the smallest scales (see iThrastarson fc Cho 



2011 



for a study of how this affects the kinetic energy spectrum of the flow). Hyperdissipation is 
applied in order to fix an artificial result of the numerical simulation and in principle should not 
have an erroneous effect on the APE generation. However, any scheme that alters the wind and/or 
temperature profiles, as this does, has the potential to change both the physical dissipation, as well 
as the efficiency with which that would be turned into reheating. 



3.2. Local theory of energetics 



In order to evaluate the importance of various possible drag mechanisms in the context of 
strongly forced atmospheres we need expressions for the local energetic s. In t his section we will 
borrow from the derivation of available potential enthalpy in iMarqueti (Il99ll ). but extended to 
explicitly track dissipation and frictional heating. 

For a steady ideal flow, without any net heating/cooling or dissipation of k inetic energ y, we 
can use Bernoulli's theorem to deflne a set of locally conserved energies (see, e.g. IVallisI |2006| ) : 



D 

m 







(11) 



where D/Dt is the material derivative, h = c^T + P/ p = CpT is the speciflc enthalpy, and $ is the 
gravitational potential energy (= gz for a constant gravitational field) Our local equations for the 



^^In Section [2] we were able to include the gravitational energy as part of the total potential energy through 
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energetics of the circulation must include these energies, the conversion terms between them, and 
local sources and sinks that can alter the energy budget locally, but should be globally conserved. 
Recognizing that not all of the enthalpy will be available for conversion to kinetic energy, we employ 
a local definition of APE as the available potential enthalpy and define it as a = /i — TrS, where 
we have subtracted off the unavailable component (UPE, u = T^s) that will globally integrate to 
agree with our previous definition of UPE (Equation [6]) . Given the definition for specific entropy 
(s = Cpln[6/9r], where the potential temperature is ^ = T[Pq/P]'^, and k = R/cp), we can expand 
out the local definition of APE: 

a = h — TrS 

= CpT - CpTr ln{T/Tr) + RTr ln{P/Pr). (12) 

Values with a subscript r are for some reference state with uniform atmospheric temperature, 
defined as: 

I _ f 1 dm 

with an implicit time-average over an appropriate interval (equivalent to the constant from 
Section [2]). The specific unavailable potential enthalpy is: 

u = TrS 

= CpTr ln{T /Tr) - RTr ln{P/Pr). (14) 





-Pr' 




P 



For an atmosphere in hydrostatic equilibrium the specific kinetic energy is dominated by the 
horizontal winds and we write: Ck = ^f^, with the understanding that v is only the horizontal 
velocity (along equipotentials). Finally, we must also track the gravitational potential energy: 
ea = gz. 

We now derive expressions for the rates of creation, conversion, and dissipation of the various 
atmospheric energies. With reference to its definition (Equation I14p . the material derivative of the 
unavailable specific enthalpy is simply: 

Du Ds (Tr\ 

m=^'m = \T)''^ ^''^ 

while from the definition of specific available enthalpy (Equation [T2|) and the thermodynamic equa- 
tion (Equation HI) we can write: 

Da _ Dh ^ Ds 
151 ~ 15i~ ''15i 



integration over the atmosphere, but since we are now defining local terms we must keep track of these components 
separately. 
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DT Tr 



" Dt 

RT 

LO h q 

P ^ 

1 /dP 



T 



dP 

— + v VP + 'u;— 
ot oz 



+ 1 



T 

± r 

Y 



1 fdP_ 
p\dt 



V- Vp] 



wg+ [1 



T 



(16) 



Here we have expanded out lo = DP/Dt = dP/dt + v ■ VP + w{dP/dz) and used the hydrostatic 
assumption: w{dP/dz) = —wpg, where w = Dz/Dt is the vertical velocity. 

For the kinetic energy we take the dot product of the horizontal velocity and the standard 
horizontal momentum equation, Dv/Dt = — (l/p)VP — f x v + Pfnc, where / = {2n sin (j))k is the 
Coriolis parameter (0 is the planetary rotation rate and (p is latitude) in the vertical direction (k), 
to get: 

DeK 



Dt 



P 



-{v-VP)-V, 



where we have defined "D as the irreversible frictional dissipation of kinetic energy (= 
The derivative of the gravitational potential energy is, by definition: 



Djec) 
Dt 



a 



D^ 
Dt 



gw 



(17) 

-V- Pfric)- 
(18) 



We explicitly split the heating into its frictional (qf) and non-frictional (mainly radiative, (?rad) 
components and combine Equations [TSlfTSl into a local energy cycle0 



Du 




'Tr' 


'm 




T 


Da 




1 - 


'm 




DeK 






Dt 




Dec 
Dt 


= wg 



(grad + Qf) 



1 



(9rad + 9/) + -(^- VP) 
P 



1 dP 

p ot 



(19) 
(20) 
(21) 
(22) 



The heating/cooling of the atmosphere ((?rad + Qf) ^cts as a source/sink of APE and UPE, with 
the local temperature structure (Tr/T) acting as an efficiency factor for generating one form of 
potential energy or the other. APE will be generated in regions of the atmosphere that are hotter- 
than-average {Tr/T < 1, so [1 — Tr/T] > 0) and heated {q > 0) or cooler-than-average {Tr/T > 1) 



^^Note that the primary difference between this energy cycle and the one presented in iMarguetl (199l|) is that we 
have exphcitly included the frictional heating as a source of potential energy. This also necessitates the inclusion of 
UPE in our set of equations. 
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and cooled {q < 0). Note that the global integral of the APE source term ([1 — Tr/T]q) is equal 

to the global source term derived previously in Equation [7] {—q*[^]*)- The APE is converted to 
kinetic energy by the term —{l/p){v- VP), which is equivalent to v ■ V{gz) if pressure is used as 
the vertical coordinate. There will be exchange between APE and gravitational potential energy 
{wg), which relates to the puffing up of a hydrostatic atmosphere when it is heated. The term 
{1/ p){dP / dt) is related to adiabatic expansion and disappears under a steady state assumption. 

If dissipated kinetic energy is returned locally as frictional heating, so that D = qf, the sum 
of Equations (|19ll22p will be D{u + a + ck + e.G)/Dt = ^radi for a steady-state atmosphere. While 
some of the processes described above may result in a local conversion of kinetic energy to heat, 
it is important to recognize that this will not always be the case. One example is magnetic drag, 
which will induce a magnetic field, and the resulting current will be dissipated non-locally. In this 
case the ohmic heating will not occur where the kinetic energy is dissipated and locally V ^ qf. 
Non-local reheating can also result when instabilities trigger gravity waves, which then propagate 



throu gh the atmosphere and deposit heat and momentum where they dissipate (IWatkins Sz Cho 



2nim . 



We can now see that the importance of returning the dissipated kinetic energy as a source of 
localized heating will depend on its strength relative to the radiative heating (gj vs. ^rad)) but also 
on its spatial correlation with the efficiency factor 1 — T,./T. Given a drag mechanism with a known 
strength and spatial distribution (and assuming localized reheating), we can use these expressions 
to estimate its effect on the energetics of the circulation. 



3.3. Drag and energetics in previously published numerical models 



Perna et al.l (l2010al ) present three numerical models that attempt to estimate the effect of 
magnetic drag on (the momentum of) hot Jupiter circulation. Although magnetic interaction 
between the winds and any planetary magnetic field could result in significant drag, the effect will 
not be localized heating since the currents sustaining the induced magnetic field will be ohmically 
dissipated in a non-local way (see lBatvgin &: Stevensonll201Cll : iPerna et al.ll2010bl . for studies of the 
related heating) . This means that we cannot use these models to estimate the amount of localized 
heating that results from magnetic drag. However, in the spirit of simplicity, we use the same drag 
strength as the magnetic drag in these models, and we assume local dissipation, in order to explore 
the effect of frictional heating on the atmospheric energetics and calibrate the level at which it may 
become important. 



We compare the models in lRauscher &: Menoul (|20ld ) and lPerna et al.l (|2010al ). which all have 
identical set-ups for a generic hot Jupiter (see the papers for details), except f or different amounts 
of Ray leigh drag applied as a term in the momentum e quation. The i nodel in iRauscher fc Menou 
(|201Cl ) has no drag applied, while the three models in iPerna et al.l (l2010al ) have drag strengths 
that vary with height (stronger higher in the atmosphere). Drag timescales range from 10^ - 10^ 
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seconds in the weakest-drag model (PMRa in Table 1), and are a factor of 10 and 100 shorter in 
the medium- and strongest-drag models (PMRb and PMRc in Table 1). For reference, advective 
timescales are on the order of 10^ seconds and radiative timescales vary from ~ 3 x 10^ to greater 
than 10^ seconds (increasing with depth). 

The local specific kinetic energy loss is calculated from the prescribed drag timescales as 
D = (l/2)w^/Tdrag- The radiative heating is treated by a Newtonian relaxation scheme, making 
it simple to determine the value of grad at each point in the modeled atmosphere. We calculate 
the mass-weighteci^ global integrals of these quantities, as well as the source terms for UPE and 
APE (from Equations [19] and [20]) , and present our results in Table 1 . Note that our numerical 
solver does not convert the lost kinetic energy into reheating (in other words, while V is nonzero 
in Equation \2T\ qj = in Equations [19] and [20]l . The values listed for the generation of UPE and 
APE are {Tr/T)q^s_d and (1 — rr./T)gi.ad) respectively, while the values listed as "missing" generation 
replace g^ad with Qf = V to determine how much potential energy should have been generated if 
we did indeed include frictional heating in the models. 

Before discussing the results, we address one issue of concern with these models, which is that 
the deepest levels (below -^10 bar) may have not cor npletely reached a s t atistic ally steady state 



by the end of the 1450 or 500 planet days simulated in lRauscher fc Menoul (120101 ) and jPerna et al 



(|2010al ). It is likely that these deep levels are still being accelerated through interaction with the 
upper atmosphere. This could complicate our energetics calculations and so we have rerun and 
extended all of the runs from those pape rs to mitigate this i s sue. We calculate energy rates at 
497 5, 5000, and 10000 planet days for the IPauscher fc MenoJ (I2OI0I ) model, and at 5000 days for 
the jPerna et al.l (j2010al ) models. 



3.3.1. Global energy rates 

The first point worth noting is that for all of these models the global integral of (/rad 7^ 0, 
meaning that the planet is not in global radiative equilibrium. There are two possible resolutions 
to this issue. The first is that the atmosphere is still being accelerated by the he ating and so has yet 



to rea ch a steady state. However, if this were the case, the value of grad for the lRauscher &: Menou 



(|2010l ) run at 10000 days should be less than at 5000 days, whereas there is not significant difference 
between them. We have also inspected the time evolution of kinetic and internal energies throughout 
the runs and we see no evidence for continued acceleration of the fiow or heating of the atmosphere. 

The second explanation, which we prefer, is that the non-zero value of grad is a measure of 
the amount of numerical dissipation in our model runs. As can be clearly seen in Equation [8] any 
amount of dissipation (physical or numerical), d, requires a non-zero average global heating rate, 
g, in order for there to be a net conversion from APE to kinetic energy to balance its loss. While 



The mass element for each point in the atmosphere is dm — pdxdydz = 



-{l/g)dxdydP. 
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Table 1. Global energetics of various runs 



RMIO^ PMRa^ PMRb'^ PMRc^ 



Radiative heating, q^ad ( x 10^^ W) 


160, 


150, 


180 


390 


430 


610 


Generation of UPE (xlO^^ W) 


-39, 


-65, 


-20 


-19 


-45 


-5.0 


Generation of APE (xlO^^ W) 


200, 


210, 


200 


410 


470 


620 


Drag dissipation, V (xlO^^ W) 









3.8 


34 


170 


Relative to grad^ 









1% 


8% 


28% 


Missing UPE gen. (xlO^^ W) 









7.3 


68 


340 


Relative to grad UPE gen.: 









38% 


150% 


6800% 


Missing APE gen. (xlO^^ W) 









-3.4 


-34 


-170 


Relative to g^ad APE gen.: 









1% 


7% 


27% 


Heat engine efficiency, r] 




7% 




13.7% 


15.7% 


20.7% 


adjusted for missing APE gen.: 








13.6% 


14.5% 


15.0% 



'^From runs identical to the one presented in iRauscher fc Menou l (I2OI0I), except 
extended out to 4975, 5000, and 10000 planet days. 



'^From runs in iPerna et al.l ( 2010al ). with increasing amounts of drag for PMRa, 
PMRb, and PMRc, associated with planetary magnetic field strengths of 3, 10, 
and 30 G, respectively. These results are for simulations run for 5000 planet days. 



Note. — All values are global totals, integrated over the atmosphere. The 
dissipation is only what has been explicitly included in the model and does not 
account for any numerical dissipation. In these models the dissipated kinetic energy 
is not re-included as localized heating, but the values we list as "missing" generation 
show what the heating would have added to the global energies. 
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physically the non-zero global heating rate will be supplied by the global frictional heating (gj in 
Equation [9]) , in our numerical models the frictional heating is not included {qf = 0) and so the 
global radiative heating, qradj must be non-zero to compensate for numerical losses (dnum)- This 
interpretation is supported by noting that q^ad is higher for the runs with drag applied, because in 
these cases the increased amount of dissipation (now both numerical and physical, dnum + ddrag) 
must be balanced by an increased amount of heating (grad) given that qf remains zero). 

In order to estimate the kinetic energy loss in our models, we can compare the global dissipation 
(assuming d = Qrad) to the input heating from incident st ellar flux. Our radia.tive fo rcing is applied 



through a Newtonian relaxation scheme (for details see iRauscher &: Menoul 120101 ) , so that there 



is no well-defined stellar flux incident on the model atmosphere. Howeve r, the values ch osen for 



the forcing are derived from a one-dimensional radiative transfer model (jiro et al.ll2005l ) for the 
hot Jupiter HD 209458b and so it is reasonable to compare our dissipation to a stellar heating 
rate appropriate for that planet Assuming an incident stellar heating rate of ~ 3 x 10^^ W, our 
drag- free model's q^ad — 1.5 x 10^^ W (= dnum) implies a numerical dissipation rate of ~ 5%. 

For the models with an increasing amount of physical drag applied (PMRa, b, c) the global 
numerical dissipation is dnum = Qvad — ^ and apparently dominates over the physical drag loss 
(D), with a mag nitude evaluated as (g^ad - T^)/{3 x 10^^) ~ 13 - 15%. Without a more detailed 
analysis it is difficult to know where this numerical loss of kinetic energy occurs in our model, 
be it evenly throughout the atmosphere or preferentially at regions of strong flow convergence or 
shear. The amount of kinetic energy lost may also be dependent on the order and magnitude 
of hyperdissipation usedEl These results imply that our numerical scheme leads to a signiflcant 
amount of kinetic energy loss for the hot Jupiter atmospheric regime, an issue that will require 
further attention outside of this paper. 

As expected, the amount of (physical) kinetic energy loss (2?) increases with the amount of 
drag applied to the atmosphere and it also becomes equal to a greater percentage of the radiative 
heating (although this is not linear with shorter drag timescales, since the increased drag will also 
work to decrease wind speeds and reduce the kinetic energy). From the formalism presented above, 
we know that it is not the ratio of Qf/qrad that matters in determining whether the energetics of 
the circulation will be altered, but the contribution that the frictional heating would make to the 
generation of APE (Equation [20]) . 

We find that in these models the spatial correlation between the atmospheric temperature 
structure and the kinetic energy dissipation is such that, in all cases, frictional heating would have 
worked to reduce APE generation, instead of feeding UPE. In these models the effect of frictional 
heating from drag on the atmospheric flow should be twofold: it will work directly to reduce wind 



^*The planet also has an internal heat flux, but that is <1% of the stellar heating. 

^^AU of the models shown here were run at T31 (spectral) resolution using hyperdissipation applied as an eighth- 
order operator on the vorticity, divergence, and temperature fields with a coefficient of 8.54 x 10*'^ m* s~^. 
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speeds and the associated heating would also work to decrease the amount of APE generated and 
available to drive the winds. 

The importance of this effect is determined by the ratio of the APE generation that would 
have come from frictional heating and the APE generation from radiative heating. We find that 
for the amount of drag applied to models PMRa, PMRb, and PMRc, the relative importance of 
APE generation from frictional heating is 1, 7, and 27%, respectively. We have demonstrated 
that — for this spatial function for the drag — the amount of drag necessary to significantly affect 
the energetics of the circulati on is someth i ng at a strength comparable to the medium-strength 



drag applied in the models of iPerna et al.l (|2010al ). We also know that the energetics would not 



have been significantly changed by the amount of drag applied to the weakest-drag model, while in 
the strongest-drag model the missing frictional heating would have significantly reduced the rate 
of APE generation. 

The energetic effect of frictional heating can also be evaluated by considering the atmosphere as 
a heat engine and calculating its efficiency. The standard Carnot efficiency is defined as rj = (Tw — 
Tc)/Tw, where T]y and Tq are the temperatures of the warm source and cold sink, respectively. We 
calculate Tw and Tc for each model as the average temperature of the regions of the atmosphere 
that are being diabatically heated and cooled, respectively, weighted by the local strength of the 
heating/cooling. This gives rj = 3, 9, 12, and 17% for RMIO, PMRa, b, and c, respectively. 
However, this is not the best estimate of the heat engine efficiency of the atmosphere. For the case 
of the Earth, the Carnot efficiency gives a value of r/ ~ 10%, but a better estimate comes from the 
actual energy input into the system: the r atio of the frictional d issipation (= APE generation) to 



incoming solar flux, which gives r] ~ 0.8% ([Peixoto &: Oortlll992l ). 



If we compare the global APE generation to incoming stellar heating (r^ 3 x 10^^ W) for our 
models we get w 7, 13.7, 15.7, and 20.7% for RMIO, PMRa, b, and c, respectively. This indicates 
that, seen as heat engines, these hot Jupiter atmospheres are at least an order of magnitude more 
efficient than the Earth's atmosphere. The efficiency increases with the amount of drag applied 
to the hot Jupiter atmosphere, as the same amount of radiative forcing has to be used more 
efficiently to generate APE, which feeds the kinetic energy, which is then dissipated at a higher 
rate. Recalculating the efficiency that would have resulted from including the "missing" APE 
generation, we get rj ~ 13.6, 14.5, and 15.0% for the PMRa, b, and c models. Since the frictional 
heating in these models is always working against APE generation, the efficiency is reduced. The 
statement that frictional heating would have reduced the APE generation by X% is equivalent to 
saying that the efficiency of the atmosphere as a heat engine would be reduced by X%. 

While this does provide an estimate of the strength of drag that may energetically affect the 
atmospheric flow — assuming local reheating from dissipation — we note that our result is dependent 
on the set-up we used, in particular the spatial form used for the drag. As we have emphasized 
throughout, it is the spatial structure of the atmosphere that dictates the flow energetics, specifically 
the correlations between temperatures, heating rates, and dissipation rates. Our results here are 
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dependent on the particular relation between the spatial function we have used for the atmospheric 
drag and the underlying atmospheric structure. Many of the drag mechanisms mentioned earlier 
will have very different spatial dependencies. The best way to estimate the energetic importance of 
any drag mechanism is to apply the Equations (jl9ti22p to a particular atmosphere with some form 
of drag at work. 



3.3.2. Spatial structure of the energetics 

We can further understand the energetics of the atmosphere by studying its spatial properties. 
This will allow us to identify the regions of the atmosphere that are the most relevant for APE 
generation, as determined by the spatial correlation between the atmospheric quantities of interest. 

Before we proceed, we need to highlight an important nuance in this analysis. Although 
Equations (jl9p and (|20p give local expressions for UPE and APE generation, note that the local 
values only have meaning in reference to the global total. The quantity Ty. is explicitly a global 
value and is an inherent part of the definitions of APE and UPE, both of which were originally 
motivated by the concept of a reference state with APE= 0. The reader should keep in mind that 
the following plots of APE and UPE generation are to be taken as an indication of which areas of 
the atmosphere (when integrated over) contribute the most to the global total. 

In Figure [T] we plot the local values of UPE and APE generation as a function of longitude 
and pressure (the model's vertical coordinate) for an equatorial slice through the atmosphere. 
(The energy rate s tend to be greates t at t he equator.) In the left panel we show the radiative 



generation for the lRauscher Menoul (120101 ) model at 10,000 days, w hile in the right pane l we plot 



the "missing" frictional generation for the strongest-drag model from iPerna et al.l (|2010al ) at 5000 
days. The pattern of radiative UPE and APE generation is similar for all of the models, although 
the amplitudes and some of the detailed structure is different. The "missing" UPE and APE 
generation in the PMRa and PMRb models is weaker and primarily concentrated in the uppermost 
layers of the atmosphere, on the night side. The longer drag times in these models mean that it 
is only the highest levels (where the applied drag is the strongest) which are able to significantly 
dissipate kinetic energy, even though most of the kinetic energy resides in the deeper layers. (See 
Figure 4 of Rauscher & Menou 2010 for a plot of kinetic energy versus depth; it has a maximum 
at ~ 2 bar.) 
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longitude [degrees east of substellar] longitude [degrees east of substellar] 



Fig. 1. — Plots of the generation of available and unavailable potential energy (APE, UPE, in watts), 
for an equatorial slice through the atmosphere (the equator dominates the energy budget). Left: 
the UPE and APE generated/lost through radiative heating/cooling in the RMIO model. Right: 
for the strongest-drag model (PMRc), the UPE and APE that would have been generated/lost 
through frictional heating, if it had been included and locally equal to the kinetic energy lost from 
drag. The colors have been scaled so that red/yellow is always positive and green/blue is always 
negative. See the text for further explanation. 



From Figure [J we can see that locally the radiative generation of UPE is much greater than 
the generation of APE, although we know from Table 1 that it will integrate to a smaller global 
rate. This is a reflection of the fact that the radiative heating on hot Jupiters is intense and a large 
amount of energy will be stored in UPE, but by its very nature this will have a small net effect on 
the atmospheric flow. 

According to the left panel of Figure [TJ most of the UPE is generated on the day side and most 
of the APE is generated on the night side (with opposite behavior for the loss). This is simple to 
explain by comparing the signs of the local efficiency factor (T^/T) and the radiative heating/cooling 
(grad)- The calculation of the reference temperature, T^, is a mass- weighted integral over the 
atmosphere (Equation [T3l) and the deepest layers, which contain the most mass, are generally 
hotter than the upper layers (see Figure 5 of Rauscher & Menou 2010 for a plot of temperature- 
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pressure profiles for different locations throughout the model). For the models under consideration 
here, ~ 1750 — 1800 and T^/T > 1 for most of the atmosphere at pressures less than 10 
bar. This means that UPE/APE will be generated where the atmosphere is heated/cooled so that 
Qrad is greater /less than zero (see Equations [T9] and [20|) . In general cool air advected from the night 
side will be heated on the day side, while flow of hot air from day to night will result in cooling 
on the night side, although we will discuss deviations from this later. The maximum UPE/APE 
generation will occur where the heating/cooling is strongest, weighted by T^/T. The amount of 
heating and cooling is a function of how much the winds are able to bring gas out of local radiative 
equilibrium and how quickly the gas can radiatively respond. These properties are themselves both 
related to the generation of APE and conversion to kinetic energy. 

The UPE and APE generation that would have resulted from localized frictional heating has 
a different spatial pattern than the radiative generation. In the right panel of Figured] we see that 
the frictional UPE generation (and APE loss) is primarily on the night side. The majority of the 
kinetic energy in the atmosphere is found on the night side and since we have adopted a horizontally 
uniform drag in these models, this results in having most of the dissipation on the night side. If 
the drag were instead stronger on the day side, for example, we might see a more even pattern of 
dissipation around the planet. Frictional heating will always be positive (the dissipation of kinetic 
energy does not lead to cooling) and so in the upper atmosphere (at P < 10 bar, where T^/T > 1) 
it will always lead to a decrease in APE. In this set-up the night side has most of the radiative 
APE generation and most of the frictional APE loss. 

As a test of how much our answers are affected by the value of the reference temperature, 
we recalculated by only integrating over the upper layers (P < 10 bar), where the atmosphere 
is being radiatively heated. This gave T,. ^ 1600 — 1650 K. We then used this value to calculate 
the global energy rates (over the entire model atmosphere) and found that our radiative APE and 
UPE generation remained unchanged, while our "missing" frictional generation was decreased by 
~ 10 — 20%. However, it is not appropriate to restrict the definition of to a partial region of the 
atmosphere. Although not directly heated, the inert layers (at P > 10 bar) contain a significant 
fraction of the total kinetic energy, which has been gained through interaction with the upper 
layers, and so must be included as part of the closed system. The real issue is then where to set to 
the bottom boundary of the model, as discussed above at the beginning of Section [3l 

We can also examine the horizontal structure of APE generation. In Figure [2] we plot the UPE 
and APE generation as a function of latitud e and longitude at the 150 mbar pressure level. The 
left panel is the radiative generation from the lRauscher fc Menoul (|20l d') model at 10, 000 days and 



the righ t panel is the "missing" frictional generation for the strongest-drag model from lPerna et al.l 



(l2010al ) at 5,000 days. The radiative generation at this level in the PMRc model is similar to 



what we see for the RMIO model, although with the amplitudes slightly decreased and the spatial 



Tr is cooler for RMIO and hotter for the PMR runs. There was a variation of ~1% in Tr between the days 
sampled for RMIO. 
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structure more smoothed out. We choose to show this particular level because it contains significant 
generation of potential energy, from both the radiative forcing and the "missing" frictional heating. 




Fig. 2. — Cylindrical projections of the generation of available and unavailable potential energy 
(APE, UPE, in watts), for a horizontal slice of the atmosphere at 150 mbar, with the substellar 
point at (0,0). Left: the UPE and APE generated/lost through radiative heating/cooling in the 
RMIO model. Right: for the strongest-drag model (PMRc), the UPE and APE that would have 
been generated/lost through frictional heating, if it had been included and locally equal to the 
kinetic energy lost. The colors have been scaled so that the red/yellow is always positive and the 
green/blue is always negative. See text for further explanation. 

As we discussed above, the generation of potential energy depends on the temperature structure 
(by Tr/T), spatially correlated with the radiative or frictional heating rate. The radiative heating 
generally works to make the day side hotter and the night side colder, with a dependence on the 
advected temperature structure, while the frictional heating depends on the flow structure and 
wind speeds. In order to more clearly understand the patterns of potential energy generation, we 
plot the temperature and velocity structures for both models at this pressure level in Figure El 

In agreement with Figure (H we see that most of the radiative UPE generation (and APE loss) 
is on the day side, while this is reversed on the night side. This is explained, as before, by the fact 
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Fig. 3. — Cylindrical projections of temperature (color, in K) and velocity (arrows, where length 
scales with magnitude) for a horizontal slice through the atmosphere at 150 mbar, with the sub- 
stellar point at (0,0). Left: for the RMIO model. Zonal (east-west) wind speeds range from -1.6 to 
4.3 km s~^ and meridional wind speeds range from -1.7 to 2 km s^^. Right: for the PMRc model. 
Zonal wind speeds range from -2.5 to 2.6 km s~^ and meridional wind speeds range from -2.8 to 
2.8 km s~^. In both models the sound speed at this level ranges from 2.1 to 2.9 km s~^. 

that the efficiency factor, T^/T, is greater than 1 at this level and so day side heating increases UPE 
while night side cooling increases APE. The deviations from this pattern are due to the detailed 
structure of the advected temperature pattern, scaled by the local efficiency factor (at this level 
Tr/T ranges between ~2 in the coldest regions and ~1.2 in the hottest regions). In general we see 
eastward advection, so that hot air is pulled from day to night across the east terminator (at 90°) 
and cold air is pulled across the west terminator. The boundary between positive and negative 
generation is not exactly at 90° longitude because the air that was heated at the substellar point 
begins to cool even before it reaches the night side. This effect is also seen at the west terminator, 
close to the poles, where the flow is westward. 

The region of reversed sign on the night side (positive UPE, negative APE generation) at 
~ ±45° latitude is because the air there is cooler than the night side equilibrium temperature and 
so is being heated, as prescribed by our Newtonian relaxation scheme. Note that night side heating 
is not unphysical. Although there is no incident stellar light, there is still flux from the deeper, 
hotter atmospheric layers, which will prevent the night side temperatures from dropping too low. 

We also find that the strongest rates of potential energy generation are not at the sub- and 
anti-stellar points ([0°, 0°] and [180°, 0°] in longitude and latitude, respectively). On the day side 
the heating is primarily of cold air advected across the west terminator (-90°) and at this level the 
relation between advective and radiative timescales is such that the maximum heating occurs at 
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~ —45°. On the night side there are two regions of coohng: where hot air has been advected from 
the day side across the east terminator, and where air has been heated by the sh ock-hke feature 



at ~ 135°. This feature is discussed in more detail in iRauscher fc Menoul (j20ld ): it is a steady 
feature that extends across many pressure levels and is associated with hot regions, consistent with 
adiabatic heating from downward flow induced by the strong horizontal convergence. Note that in 
Figure [l] there is a region of night side APE loss associated with this feature. 

While the potential energy generation from both the radiative and frictional heating will be 
weighted by the local temperature through the efficiency factor T^/T, the radiative heating is more 
strongly dependent on the temperature structure, while the frictional heating is dependent on the 
pattern of wind speeds. The "missing" potential energy generation shown in the right panel of 
Figure [2] follows the same structure as the kinetic energy at this level. We can see this from the 
flow pattern in Figure[3l where the strongest winds are on the night side, maximized in an equatorial 
jet and with significant flow in curving branches at high latitude. 



4. A Special Case: Drag in Eccentric Planet Atmospheres 

Not all of the strongly irradiated close-in gas giants are on circular orbits; there are many known 
with s ignificant non-zero e ccentricities, the most extreme being HD 80606b, with an eccentricity of 



0.934 (|Moutou et al.ll2009l ). These eccentric planets may be subject to the same drag mechanisms 
that have been proposed for circular hot Jupiters, but with the added complication that the level of 
stellar irradiation will change throughout the orbit. During periastron passage there will be intense 
stellar heating that will drive up the potential energy of the atmosphere, some of which will be 
converted into kinetic energy and then dissipated. The temperature and velocity structures of the 
atmosphere will vary throughout the orbit. The radiative APE generation should closely follow the 
temperature evolution of the atmosphere, while the frictional dissipation will follow the evolution 
of the kinetic energy. At any point in the orbit there may be a mismatch between the rates of 
APE generation, conversion to kinetic energy, and dissipation, leading to a temporary build up of 
APE or kinetic energy. This delay in energy transfer is seen in the three-dimensi onal model of the 



eccentric hot Neptune GJ436b (a = 0.02872 AU and e = O.lS. lLewis et al.ll20ld ). where they find 
that the peak in wind speeds occurs 4-8 hours after periastron passage. The evolving dissipation 
may then lead to a secondary period of APE generation from the frictional heating and, depending 
on the delayed rates of energy transfer, this could be late enough after periastron passage that the 
frictional heating might dominate over the stellar heating. As we can see, steady state assumptions 
and implicit time-averages are no longer valid for eccentric planet atmospheres and a more detailed 
analysis of the energetics will require further study. 
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5. Summary and Conclusions 

We have presented a discussion of the standard formahsm for understanding atmospheric 
energetics, exphcitly accounting for frictional heating from the dissipation of kinetic energy. The 
potential energy of the atmosphere is divided into the component that can be converted into 
kinetic energy (the "available" potential energy, APE) and the component associated with uniform 
changes in temperature (the "unavailable" potential energy, UPE). Uniform frictional heating will 
feed the UPE, while differential heating (from any source) will feed APE. Frictional heating can 
be safely neglected from the energetics of the atmosphere (to some level of precision) if it is: 1) 
predominantly uniform, or 2) its differential heating is much less than the radiative differential 
heating and it contributes minimally to the generation of APE. 

We give local expressions for energy generation, conversion, and dissipation in Equations (jlQl 
[22]) . APE is generated in regions with positive correlations between temperature and heating 
(when hotter-than-average regions are heated and colder-than-average regions are cooled) . In order 
to determine if frictional heating matters for the energetics of an atmosphere, we must identify 
possible drag mechanisms and estimate their strengths and spatial variations. We can then take 
global integrals of Equations P9II22|) to determine if the APE generation from frictional heating 
is a significant fraction of the radiative APE generation. Given the many unknowns related to 
modeling hot Jupiter atmospheres, one may consider that frictional APE generation is no longer 
negligible once it exceeds ~10% of the radiative APE generation. This is equivalent to saying that 
the heat engine efficiency of the atmosphere will be altered by more than 10%. 

As an exercise and test of this formalism, we calculated APE and UPE generation rates for 
numerical models that include drag as a kinetic energ y sink but do not return th e en ergy as heat 



We use d the previously published numerical models bv lRauscher &: Menoul (j2010|) and lPerna et al 



(|2010al ). which have identical set-ups but a range of representative drag strengths. We found that 
in these models the spatial correlations between temperature and heating rates are such that in 
all cases the frictional heating would have worked to decrease the global generation of APE (or, 
equivalently, the atmosphere's heat engine efficiency), at a factor of of 1, 7, and 27% for the models 
with weak, medium, and strong drag, respectively. This provides an estimate of the drag strength at 
which frictional heating can significantly alter the energetics of the circulation, with the caveat that 
this result is strongly dependent on the spatial form we chose for the applied drag. In these models 
drag timescales were horizontally constant; drag mechanisms that have greater spatial variation 
can have more effect on the energetics at weaker drag strengths. 

Through this analysis we were also able to estimate the rate of numerical dissipation in these 
models. Any atmospheric dissipation of kinetic energy (numerical or physical) must be balanced 
by a non-zero global heating rate, in order to generate the APE that is converted to kinetic energy 
(Equation [8]) . In these models the only source of heating is radiative and so this heating must be 
non-zero, although it means that the atmosphere is not in global radiative equilibrium. Comparing 
the net radiative heating rate to a representative value for the incident stellar flux, we estimate 
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that the rate of numerical kinetic energy loss from these models is substantial, at 5 — 15%. 

By calculating local APE generation throughout these models we were able to identify the 
regions that contributed the most to the global energy rates. The source term for APE in Equation 
(I20p is a function of the local heating rate and an efficiency factor based on the local temperature 
structure, 1 — T^/T, where T^- is a mass- weighted harmonic mean of temperature throughout the 
entire atmosphere (Equation [l3]). The deeper layers of the atmosphere, which contain most of 
the mass, are generally hotter than the upper atmosphere. In these models T^/T > 1 for pressures 
less than ~10 bar, the same region in which all of the stellar radiative heating occurs. This means 
that at these pressures APE is generated radiatively where the atmosphere is cooled (mostly on 
the night side) and is lost where the atmosphere is heated (mostly on the day side). 

In these models most of the kinetic energy resides on the night side and — with our choice of a 
horizontally constant drag — this means that most of the dissipation occurs on the night side. Since 
frictional dissipation always produces heating (never cooling) , in this model frictional heating always 
decreases the APE at pressures less than 10 bar. Under our assumption that any lost kinetic energy 
is locally returned as heating, this means that the APE loss due to friction is spatially correlated 
with the regions of strong kinetic energy and strong radiative APE generation. 

Finally, we discuss how the energetics will change for planets on eccentric orbits. Since the 
temperatures and wind structure of the atmosphere should vary throughout the planet's orbit, 
we expect evolving rates of APE production, conversion, and dissipation of kinetic energy. A 
mismatch, or time delay, between these energy transfer rates would result in a build up of APE or 
kinetic energy, and following periastron passage there could be a secondary period of heating from 
frictional dissipation. 

Through this analysis we have shown how frictional heating would alter the generation of 
available potential energy and have argued that if this effect is substantial, it could change the 
nature of the atmospheric circulation, since the APE is converted into kinetic energy (in the form 
of winds). However, it is difficult to predict what effect we would see in the temperature and 
wind structure of the atmosphere. In order to test how big of a change would result from the 27% 
decrease in APE generation for the PMRc model, for example, we would need to run a model that 
explicitly included the localized heating from frictional dissipation. While beyond the scope of this 
paper, this will be an interesting avenue to pursue in future work. 

In conclusion, we suggest that performing an energetic analysis of atmospheric models for 
strongly-irradiated planets can be beneficial in several ways. Aside from providing another tool 
with which to understand these exotic atmospheres, the diagnostics presented here can be used 
to evaluate whether particular heat sources (such as frictional dissipation) can significantly affect 
the atmospheric circulation and should be explicitly included in numerical models. Additionally, 
this analysis provides an estimate of the amount of numerical kinetic energy loss in models. While 
there should be some amount of physical dissipation in an atmosphere, we need to carefully consider 
whether the magnitude and spatial properties of the numerical dissipation introduce any erroneous 
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effects, and further analysis is warranted. 

We thank Adam Showman and Jeremy Goodman for valuable comments that helped improve 
the quality of this paper. We thank the anonymous referee for useful feedback. This research 
began while the authors were in residence at the Kavli Institute for Theoretical Physics, generously 
supported by the National Science Foundation under Grant No. PHY05-51164. This work was 
performed in part under contract with the California Institute of Technology (Caltech) funded by 
NASA through the Sagan Fellowship Program. KM recognizes support from NASA under grant 
no. PATM NNX11AD65G. 

REFERENCES 

Batygin, K., & Stevenson, D. J. 2010, ApJ, 714, L238 

Dobbs-Dixon, I., Gumming, A., & Lin, D. N. G. 2010, ApJ, 710, 1395 

Dobbs-Dixon, I., & Lin, D. N. G. 2008, ApJ, 673, 513 

Button, J. A. 1973, Tellus, 25, 89 

Goodman, J. 2009, ApJ, 693, 1645 

Iro, N., Bczard, B., & Guillot, T. 2005, A&A, 436, 719 

Fortney, J. J., & Nettelmann, N. 2010, Space Sci. Rev., 152, 423 

Kalnay, E. 2002, Atmospheric Modeling, Data Assimilation and Predictability, by Eugenia Kalnay, 
pp. 364. ISBN 0521791790. Cambridge, UK: Cambridge University Press, December 2002 

Landau, L. D., & Lifshitz, E. M. 1959, Fluid Mechanics, Oxford: Pergamon Press, 1959 

Lewis, N. K., Showman, A. P., Fortney, J. J., Marley, M. S., Freedman, R. S., & Lodders, K. 2010, 
ApJ, 720, 344 

Li, J., & Goodman, J. 2010, ApJ, 725, 1146 

Lorenz, E. N. 1955, Tellus, 7, 157 

Marquet, P. 1991, Quarterly Journal of the Royal Meteorological Society, 117, 449 
Menou, K., k Rauscher, E. 2009, ApJ, 700, 887 
Moutou, C., et al. 2009, A&A, 498, L5 

Pearce, R. P. 1978, Quarterly Journal of the Royal Meteorological Society, 104, 737 



-25- 



Peixoto, J. P., & Oort, A. H. 1992, Physics of Climate, by Jose P. Peixoto & Abraham H. Oort, 
New York: American Institute of Physics (AIP), 1992 

Perna, R., Menou, K., & Rauscher, E. 2010, ApJ, 719, 1421 

Perna, R., Menou, K., & Rauscher, E. 2010, ApJ, 724, 313 

Rauscher, E., & Menou, K. 2010, ApJ, 714, 1334 

Showman, A. P., Cho, J. Y.-K., & Menou, K. 2010, in Exoplanets, ed. S. Seager, 471 

Showman, A. P., Cooper, C. S., Fortney, J. J., & Marlcy, M. S. 2008, ApJ, 682, 559 

Showman, A. P., Fortney, J. J., Lian, Y., Marley, M. S., Preedman, R. S., Knutson, H. A., k, 
Charbonneau, D. 2009, ApJ, 699, 564 

Showman, A. P., & Guillot, T. 2002, A&A, 385, 166 

Thrastarson, H. T., & Cho, J. Y-K. 2011, ApJ, 729, 117 

ValHs, G. K. 2006, Atmospheric and Oceanic Fhiid Dynamics, by Geoffrey K. VaUis, pp. 770. Cam- 
bridge University Press, November 2006. ISBN-10: 0521849691. ISBN-13: 9780521849692 

Watkins, C, & Cho, J. Y.-K. 2010, ApJ, 714, 904 



This preprint was prepared with the AAS IM^^X macros v5.2. 



